Spin glasses in the limit of an infinite number of spin components. 
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We consider the spin glass model in which the number of spin components, m, is infinite. In 
the formulation of the problem appropriate for numerical calculations proposed by several authors, 
we show that the order parameter defined by the long-distance limit of the correlation functions is 
actually zero and there is only "quasi long range order" below the transition temperature. We also 
show that the spin glass transition temperature is zero in three dimensions. 

PACS numbers: 



I. INTRODUCTION 

It is of interest to study a spin glass model in which 
the number of components of spin components m is in- 
finite, because it provides some simplifications compared 
with Ising (m = 1) or Heisenberg (to = 3) models. For 
example, in mean field theory (i.e. for the infinite range 
model) there is no "replica symmetry breaking" i so the 
ordered state is characterized by a single order parameter 
g, rather than by an infinite number of order parameters 
(encapsulated in a function q{x)) which are needed^ for 
finite-TO. 

There has recently been renewed interest in the 
TO = oo model, and the interesting result emerged from 
these studies that the effective number of spin compo- 
nents depends on the system size N and is only really 
infinite in the thermodynamic limit. One motivation for 
the present study is to investigate some consequences of 
this result. 

Further motivation for our present study comes from 
earlier work by two of us^ which argued that an isotropic 
vector spin glass (to > 2), as well as an Ising spin glass, 
has a finite spin glass transition temperature Tsg in three 
dimensions. The results of Ref. Q also indicate that 
Tsg is very low compared with the mean field transition 
temperature, T^'q^ , and decreases with increasing m, see 
Table n The results in Table |I| suggest that Tsa/T^J 
may be zero in the m = oo limit in three dimensions, and 
we investigate this possibility here. 

In this paper, we study the to = oo SG model, both the 
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TABLE I; Estimates of the spin glass transition temperature, 
relative to the mean field value, Tsg^ = y/z/m, see Eq. l(2j, 
for different values of m for the three-dimensional simple cubic 
lattice (2 = 6). The factor of l/m in T^r f ears because 
the spins were normalized to unity in Refs. |5|6| |7||, rather than 



infinite range version and the short-range model in three 
and two dimensions. We find that we need to carefully 
specify the order in which the limits m — > 00 and the 
thermodynamic limit N ^ 00 are taken. In Ref. the 

— > 00 limit is taken first (since a saddlepoint calcula- 
tion is performed) and the to ^ 00 limit is taken at the 
end. However, in the formulation of the m = 00 problem 
which has been proposed for numerical implementation 
in finite dimensions'^'tttR and which we use here, the limit 
TO —> 00 is taken first for a lattice of finite size. In the 
latter case, we find that for T < Tsg the spin glass cor- 
relations decay with a power of the distance r and tend 
to zero for r 00, so the order parameter, defined in 
terms of the long-distance limit of the correlation func- 
tion, is actually zero. Nonetheless, there can still be a 
transition at Tsg separating a high temperature phase 
where the correlations decay exponentially, from the the 
low temperature phase where they decay with a power 
law. By contrast, if one takes N 00 first with to finite, 
the power law decay eventually changes to a constant at 
large r and so a non-zero spin glass order parameter can 
be defined, as in Ref. U 

We give phenomenological arguments for these conclu- 
sions and back them up (for the case where to —> 00 is 
taken first) by numerical results at zero temperature. We 
also find, from numerical results at finite temperature, 
that Tsg /Tsg' = in three dimensions for to = 00, 
consistent with the trend of the results in Table HJ 

In Sec.niwe discuss the model and the methods used to 
study it numerically. In Sec. lIIII we describe our results at 
T = for both short-range and the infinite-range model, 
while in Sec. II VI we describe finite temperature results for 
short-range models. Our conclusions are summarized in 
Sec. El 



II. MODEL AND METHOD 

We take the Edwards- AndersoniSi Hamiltonian 



(1) 



to m}^'^ as here. 
finite for M —> o 



For the model used in this paper, T^q is 



where the spins Si {i — 1,---,N) are classical vectors 
with TO components and normalized to length m}^'^, i.e. 



2 



Sf = m. As we shall see, this normalization is necessary 
to get a finite transition temperature in the mean field 
limit. The Jy are independent random variables with 
a Gaussian distribution with zero mean. We consider 
both the infinite range model and short-range models 
with nearest-neighbor interactions in two and three di- 
mensions. For the infinite range model, the standard 
deviation is taken to be — 1 while for the short- 

range models the standard deviation is set to be unity. 
According to the mean field approximation, the spin glass 
transition temperature is 
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(2) 



where [• • -Jav indicates an average over the disorder. 
Hence, for the infinite range model, (where mean field 
theory is exact) Eq. © gives Tsg = Tgcf = 1' "^^^^^ 
the short-range case it gives Tgjf — y^, where z is the 
number of nearest neighbors (4 for the square lattice and 
6 for the simple cubic lattice). 

As shown in other work2iiiSi£j the problem can be sim- 
plified for m ^ oo. The spin-spin correlation function. 



ij — —{Si ■ Sj), 



is given by 



T-^a, = (A^'),^., where 



(3) 



(4) 
(5) 



and the Hi have to be determined self consistently to 
enforce (on average) the length constraint on the spins, 



Ciq 



1 



(6) 



Angular brackets, (•••), refer to a thermal average for a 
given set of disorder. Eq. © with i = 1, • • • , TV represents 
N equations which have to be solved for the N unknowns 
Hi. In Sec. II VI we will solve these equations numerically 
for a range of sizes at finite temperature. We emphasize 
that in Eqs. ®~® the limit m ^ oo has been taken 
with TV finite. This is the opposite order of limits from 
that in the analytical work of Ref. Q where N ^ oo was 
taken before m — > oo. As we shall see, the results from 
the two orders of limits are different. 



Eqs. are not well defined at T = 0. However, 

Aspelmeier and Moore* pointed out that one can solve 
the TO = oo problem directly at T = 0, using the follow- 
ing method. At zero temperature there are no thermal 
fluctuations so each spin lies parallel to its local field, i.e. 



H, 



(7) 
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FIG. 1: A plot of the effective number of spin components as 
a function of system size A*'. For small TV, mo ~ TV, but once 
mo hits the actual number of spin components m, it sticks at 
m as TV is further increased. 



where m^^^Hi is the magnitude of the local field on site i. 
Remarkably, it was shown by Hastings^ that these local 
fields are precisely the zero temperature limit of the Hi 
in Eq. Hastings^ also showed that the number of 

independent spin components which are non-zero in the 
ground state, which we call toq, cannot be arbitrarily 
large, but satisfies the bound 



Too < 



(8) 



This means that one can always perform a global rota- 
tion of the spins such that only toq components have a 
non-zero expectation value and the remaining to — toq 
components vanish. Thus one can think of toq as the 
effective number of spin components. If to is finite, then, 
at some value of TV, toq would equal the actual number of 
spin components to. At this point, all spin components 
are used so too "sticks" at the value to as TV is further 
increased, see Fig. ^ 

More generally we can write Eq. © as 



Too ~ TV^, (too < to) 



(9) 



and the bound in Eq. ((HJ gives fj, < 1/2. Later, we will 
determine /i numerically for several models. For Eqs. Q- 
® to be valid we need to > too which corresponds to the 
curved part of the line in Fig.^ As discussed above, this 
corresponds to taking the limit to — > oo first, followed by 
the limit TV — > oo. Since too increases with TV one needs 
larger values of to for larger lattice sizes. This will be 
important in what follows. 

We therefore see that we can numerically solve the to — 
oo problem at T = on a finite lattice by taking a number 
of spin components which is finite but greater than too, 
and solving Eqs. JT)). To do this we cycle through the 



3 



lattice, and at site i, say, we calculate Hi from 



m 



1/2 



(10) 



We then set S,; to the value given by Eq. {Tj) so it lies 
parallel to its instantaneous local field. This is repeated 
for each site i, and then the whole procedure iterated it- 
erated to convergence. Although spin glasses with finite- 
m have many solutions of Eqs. iQ, it turns out that for 
m = oo (in practice this means m > mo) there is a unique 
stable solutionii, so the numerical solution of Eqs. (0 is 
straightforward. We will discuss our numerical results 
at T = using Eqs. JT)) in Sec. IIIII and here we simply 
note that we do indeed find a unique solution of these 
equations. 

Next we consider the order parameter in spin glasses 
for m = oo. In the absence of a symmetry breaking field, 
one defines the long range order parameter, q, by the 
behavior of the spin-spin correlation function [Cj%j]av at 
large distances, i.e. 



q2= hm [Cf,]a 



(short range), 



(11) 



where Rij = jR^ — Rj|. For the infinite-range model, any 
distinct pair of sites will do, and so 



q 



(« ^ ]) (infinite range). (12) 



We now give phenomenological arguments, which will 
be supported by numerical data in Sec. IIIII that q ob- 
tained from Eqs. and (|12|l . in which Cy is deter- 
mined by Eqs. Q-JSl), is actually zero for m = oo, 
and that, at best, spin correlations have only "quasi- long 
range order". For the short range case, this means that 
[Ci'^]av decays with a power of the distance Rij, while 
for the infinite range model the correlation function in 
Eq. H12|) tends to zero with a power of A^. 

To see why this is the case, we take T = and consider 
first the infinite-range model. For a given TV, the spins 
"splay out" in mo ^ A^ directions. We expect the spins 
to point, on average, roughly equally in all directions in 
this mo-dimensional space. Now Cij in Eq. ^ is equal to 
cos 9ij where 6ij is the angle between Si and Sj. We take 
the square and average equally over all directions. To do 
the average, take a coordinate system with the polar axis 
along Si , so 6*^ = 9j the polar angle of Sj . Then we have 



(COS 



1 E'^liisl) 1 



mo 



= N-f , (13) 

mo 



where we used the result that the average is roughly the 
same for all the too spin components. Since fi will turn 
out to be non-zero it follows that the order parameter 
tends to zero with a power of the size of the system. 
The same will be true at temperatures T < Tsg, while 



above Tsg the order parameter as defined here will vanish 
faster, as 1/A^. 

How can we reconcile this vanishing order parameter 
with earlier results^ that the order parameter is non zero 
below Tsg — 1; and in particular is unity at T = 0. The 
difference comes in part because q^ in Ref . 1] , which we 
call QajkT' times our q^, and so 



9AJKT = mq 



o 9 TO , , 

.2 2 _ ^ T= 0). 

mo 



(14) 



The other difference is that Ref. l] performs the limit 
N oo first, which corresponds to being on the horizon- 
tal part of the line in Fig. ^ so toq = m. From Eq. H14(l 
we then get (/aikt ~ const, at T = in agreement with 
Ref. 

Going back to the calculation of Cij , if one sums [C|^ ]av 
for the infinite range model over all pairs of sites we find 
that the spin glass susceptibility xsG at T = is given 
by 



XSG 



= ^ Et^^lav = 1 + (iV - 1)9' ^ Nq' ^ N'^ 



(15) 

Turning now to the short-range case, we expect that 
XSG ~ N^^'^ will still be true, which implies that corre- 
lations decay with a power of distance. Assuming that 
[Cfj]e,v ^/^ij some exponent y, then integrating 
over all r up to r = i (where A^ — U^) and requiring 
that the result goes as N^~^, gives y — d fi, i.e. 



1 



.d lA 



(16) 



Such power law decay is often called "quasi long range 
order". We expect that Eq. (|16|) will be true quite gen- 
erally at T = and everywhere below Tsg if Tsg > 0. 
Note that this implies that q = according to Eq. Hll|l . 
Above Tsg, [Cfj]iiv will decay to zero exponentially with 
distance. 

If TO is large but finite, then [CJ^^■]av will saturate when 
Rij is sufficiently large that all the spin components are 
used. This happens when [CJ^^]av ~ 1/to, i.e. for Rij > 
^i/dfi^ In this case, (/ajkt = "^9'^ '^i^^ be finite according 
to Eq. ((TTl). 



In Sees, nil XI and ITll Bl we will provide numerical sup- 
port for Eq. 11511 for the infinite-range and short-range 
cases respectively. 



III. RESULTS AT ZERO TEMPERATURE 
A. Infinite Range Model 

We consider a range of lattice sizes up to A = 2048 and 
for each size the number of samples is shown in Table ITU 
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N 




32 


1000 


64 


1000 


128 


1000 


256 


1000 


512 


1000 


1024 


777 


2048 


302 



TABLE II: Number of samples used in the T = studies of 
the infinite-range model. 
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FIG. 2: The average number of non-zero spin components in 
the ground state, mo, as a function of A'^ for the infinite-range 
model. We see that mo increases like A'' with fj. close to 2/5 
as expected. 

The average number of non-zero spin components in 
the ground state is given by Eq. for which it has 
been shown thati'** 

fi = 2/5 , (infinite range.) (17) 

exactly. This result has been confirmed numerically^. 
Our results for /i are shown in Fig. |21and indeed give ^ 
close to 2/5. The small deviation is presumably due to 
corrections to scaling. 

We also calculated at T = from Eq. ifT^ . In 
Eq. (PI the thermal average, (•••), is unnecessary, and 
the spin directions are determined by solving Eqs. Q) 
and pU|l . The results for arc shown in Fig. O showing 
that it vanishes with exponent — /i as a function of N, as 
expected from Eq. ifT^ . 




102 IQ3 
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FIG. 3: The square of the order parameter at T = for 
infinite range model. As expected it decreases like A^^ with 
H = 2/5. 



B 6 



2 





I 1 I 1 1 1 1 1 1 1 1 1 1 1 M 1 

slope = /J, = 0.33 y« 


I 1 1 1 




y 3d 






X T = 






1 1 1 1 1 1 1 1 1 1 1 1 


I 1 1 1 



10 102 IQ3 

N 



FIG. 4; The average number of non-zero spin components 
in the ground state, mo, as a function of N for the short- 
range model in d = 3. We see that mo increases like A*" with 
11 ~ 0.33. 
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TABLE III: Number of samples used in the calculations for 
the short-range model in three-dimensions. 
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FIG. 5: The spin glass susceptibility for the short-range 
model in d = 3 for different system sizes. As expected it 
varies as A^^"*", where fi ~ 0.33 was also found in Fig. 0] 
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TABLE IV: Number of samples used in the calculations for 
the short-range model in two-dimensions. 
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B. Short Range models 



FIG. 6: The average number of non-zero spin components 
in the ground state, mo, as a function of A' for the short- 
range model in d = 2. We see that mo increases like A''' with 
fj, ~ 0.29. 



First of all we describe our results for three dimensions. 
The number of samples is shown in Table Hill 

Our results for fj, are shown in Fig. ^ indicating that 
fi ~ 0.33, definitely different from the infinite range result 
of 2/5. The results for xsG as a function of N are shown 
in Fig. |31 We see that xsG grows with an exponent 1 — /i 
with the same value of /x as in Fig. 0] We therefore find 
that ~ 1.0, and so, from Eq. H16|) . the spin glass 



correlations decay as 

[C^^l-^, (d = 3, r = 0). (18) 

(It is of course possible that power of Rij may not be 
exactly —1.) 

Next we describe our results for two dimensions. The 
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FIG. 7: The spin glass susceptibility for the short-range 
model in d = 2 for different system sizes. As expected it 
varies as hS^^^\ where p ~ 0.29 was also found in Fig.|S| 



number of samples used is shown in Table llVl Our results 
for /i are shown in Fig. El and give \x ~ 0.29. The data 
for xsG is shown in Fig. [3 We see that XSG increases as 
7V^~^ with the same [i as determined from Fig. |H1 We 
therefore find that d/i ~ 0.58, and so, from Eq. (|16() . the 
spin glass correlations decay as 



nO.58 ' 



(d = 2, T = 0) , 



(19) 



IV. 



RESULTS FOR SHORT RANGE MODELS 
AT FINITE TEMPERATURE 



We have determined finite temperature properties by 
solving Eqs. Q-® self-consistently using the Newton- 
Raphson method. We start at high temperature, T = Ti 
say, and take our initial guess to be iJ^ = 1//3 which 
is the solution obtained perturbatively to first order in 
1/T. We then solve the equations at successively lower 
temperatures, Ti > T2 > > T-^ ■ ■ and obtain the ini- 
tial guess for the Hi at temperature Ti^i by integrating 
the equations"' 
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FIG. 8: The spin glass susceptibility as a function of temper- 
ature in three dimensions. The vertical axis has been divided 
by j^d(i-M)^ in which we took /i = 1/3 in order to collapse the 
data at T = according to the data in Figs. |1] and |3 
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FIG. 9: A scaling plot of the spin glass susceptibility in Fig.|Hl 
assuming a zero temperature transition. 



from ft toft+i (/3 = l/T). 

Results for xsG in d = 3 are shown in Fig. |H1 in which 
we scaled the vertical axis by L'^^^-t^) L^) so the data 
collapses at T = 0. If we assume a zero temperature 
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FIG. 10: The main figure is a plot of /L against T in three 
dimensions. The inset shows ^l/L at T = as a function of 
L. The dashed line is a guide to the eye. 
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FIG. 11: A scaling plot of the data for XSG in two dimensions, 
assuming a zero temperature transition. In the vertical axis, 
XSG is divided by ~ L^-^^ so that the data collapses 

at T = 0. 
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transition, the data should lit the finite-size scaling form 

XSG = L'^^^-^^X [l^'^T^ . (22) 

where X[x) —>■ const, for x 0, and the power law 
prefactor in front of the scaling function X{x) then gets 
the T = hmit correct. Figure El shows an appropriate 
scaling plot with i/ = 1.23. Apart from the smallest size, 
L = 4, the data clearly collapses well. By considering 
different values of v we estimate 



v = 1.23 ±0.13 



(d = 3) , 



(23) 



This result can be compared with that of Morris et ali^ 
who quote v = 1.01 ± 0.02. Since our results cover a 
larger range of sizes and have better statistics, we feel 
that the error bars of Morris et al. are too optimistic. 
Assuming this, our result is consistent with theirs. 



intersect at Tsg and also splay out below Tsg- To deter- 
mine we Fourier transform [CJ^^■]av to get xsG(k) and 
then use^*i^ 
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/ Xsg(O) 

2sin(fcmi„/2) V 

XSG(kmin) 



- 1 



1/2 



(25) 



where 'kmin = (27r/L)(l, 0, 0) is the smallest non-zero 
wavevector on the lattice. 

The results are shown in the main part of Fig. 1101 
The data don't intersect at any temperature, but seem 
to be approaching an intersection at T = for the larger 
sizes. To test out this possibility, we have computed the 
correlation length directly at T = 0, from the solution 
of Eqs. (TJ) and (|10|l . where we can study larger sizes 
than in the finite-T formulation of Eqs. Q-lO. The 
data is shown in the inset to Fig. ^| It indicates, fairly 
convincingly, that ^l/L approaches a constant for L — > 
CXI at T = 0, and hence that there is a transition at T = 0. 



We should, however, also test to see if the data can 
be fitted with a finite value for Tsg- To do this, it is 
convenient to analyze the correlation length of the finite 
system, , and plot the dimensionless ratio ^l/L which 
has the expected scaling form^ii^ 

^^f(^L'/''{T-Tsg)) (24) 

without any unknown power of L multiplying the scal- 
ing function F. Hence the data for different sizes should 



In d = 2 it is well established that Tsg = even for 
the Ising case. A scaling plot for xsG for m = oo in 
d = 2, corresponding to Eq. 122|) . is shown in Fig. ^] 
with 1/ — 0.72, which gives the best data collapse for 
larger sizes, and d{l — /i) = 1.42 which is obtained from 
the T = results in Sec. IIIII Again the data scales well. 
Overall we estimate 

v = 0.72 ± 0.05 (d = 2, from xsg) • (26) 




FIG. 12; Data for jL as a function of T in two dimensions. 
Clearly the data for larger sizes is merging at T = indicating 
a transition at Tsg = 0. The inset shows data for ^l/L at 
T = confirming that the data becomes independent of size 
at T = 0. The dashed line is a guide to the eye. 

This is consistent with the results in Morris et ali^i who 
quote V = 0.65 ± 0.02. 



We have also computed the correlation length in 
two dimensions, and show the data in Fig. (|12|l . The 
curves become independent of size, for large L, at T — 0, 
confirming that Tsg = 0. A scaling plot of the data for 
the largest sizes {L > 24) in Fig. ^| has the best data 
collapse with ly — 0.65 and altogether we estimate 

= 0.65 ±0.05 (rf = 2, from^L/i), (27) 

which is consistent with our estimate from xsG in 
Eq. 1261) , and with the result of Morris et ali^. 



V. CONCLUSIONS 

We have considered the spin glass in the limit where 
the spins have an infinite number of components. In 
the formulation of this problem appropriate for numer- 
ical calculations^iiiSiS, where the limit m — > oo is taken 
with N finite, we find that the order parameter, de- 
fined in terms of correlation functions in zero (symmetry- 
breaking) field, vanishes. Instead, below Tsg-, there is 
only "quasi-long range order" in which the correlations 



FIG. 13: A scahng plot of the data from the largest sizes for 
^l/L in two dimensions assuming Tsg = 0. 



decay to zero with a power of distance. Whereas we 
define the order parameter in terms of the the long dis- 
tance limit of the correlation functions, Aspelmeier and 
Moor^ define a local order parameter in terms of the 
contribution to the constraint in Eq. @ that comes from 
the eigenmodes with zero eigenvalue of the matrix Aij. 
They argue their order parameter is related to the sus- 
ceptibility in the presence of a small field h, where the 
limit A'^ — > CX3 is taken before the limit h —* in order to 
break the symmetry. From numerics on the infinite-range 
model, Aspelmeier and Moore claim that their order pa- 
rameter agrees with that of Almeida et alA. 

However, in a sensible physical model, any reasonable 
definition of the order parameter should give the same 
answer. In particular, one should be able to obtain the 
square of the order parameter from the long distance 
limit of the correlation function (off-diagonal long range 
order) in zero field, and get the same answer as the local 
expectation value of the spin in the presence of a small 
symmetry breaking field. This does not appear to be the 
case for the m — oo model if the limit m — > oo is taken 
before N ^ oo. 

On the other hand, if the thermodynamic limit, N — > 
oo, is taken with m large but finite, then the correlations 
saturate at a value of order 1/m at large distance, and so 
a finite spin glass order parameter can be defined from 
the long distance limit of the correlation functions. This 
seems to agree with that found in the analytical work of 
Ref. 0, and is presumably the same as the local order 
parameter in a symmetry breaking field. Hence, there 
seems to be no inconsistency if the limit A — > oo is taken 
first. 
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We have also studied the m = oo model in three di- 
mensions, finding the transition to be at zero tempera- 
ture, in contrast to the situation iov^^ m = 1,2 and 
3. We suspect that Tsg = only in the m = oo limit, 
rather than for all m less than some (non-zero) critical 
value nic, since spin glasses with m = oo seem to have 
unique features. We have already mentioned that there 
is only quasi long-range order below Tsg in this case, 
in contrast to finite-m. Another example is that Green 
et al^ find the upper critical dimension, above which 
the critical exponents are mean field like, to be du = 8, 
whereas for finite m one has (i„ = 6. Our result that 
Tsg = for m = oo in d = 3 is consistent with the claim 
of Viana™ that that the lower critical dimension (below 
which Tsg — 0) is also di = 8, but currently we cannot 
say anything specific about dimensions above 3. 

We find, not surprisingly, that Tsg — also in two 
dimensions. Our results for the correlations length expo- 
nent at the T — transition in d = 2 and 3 are consistent 
with those of Morris et al^. 

Finally, we note that Aspelmeier and Moor^ have pro- 



posed that the m = oo model is a better starting point 
for describing Ising or Heisenberg spin glasses in finite di- 
mensions, than the Ising model. We have argued in this 
paper that the spin glass with m strictly infinite is not a 
sensible model, but one rather needs to consider m large 
but finite. Hence, the m = oo formulation proposed by 
Aspelmeier and Moore'* and othersiSiS would need to be 
extended to a 1/m expansion and evaluated, at the very 
least, to order l/m. More probably an infinite resumma- 
tion would be needed (M. A. Moore, private communi- 
cation) to obtain sensible results in the spin glass phase, 
but this may be feasible. 
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